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We describe the results of a systematic high-statistics 
Monte-Carlo study of finite-size effects at the phase transition 
of compact U(l) lattice gauge theory with Wilson action on 
a hypercubic lattice with periodic boundary conditions. We 
find unambiguously that the critical exponent v is lattice-size 
dependent for volumes ranging from 4 4 to 12 4 . Asymptotic 
scaling formulas yield values decreasing from v{L > 4) « 0.33 
to v{L > 9) ~ 0.29. Our statistics are sufficient to allow the 
study of different phenomenological scenarios for the correc- 
tions to asymptotic scaling. We find evidence that corrections 
to a first-order transition with v = 0.25 provide the most ac- 
curate description of the data. However the corrections do 
not follow always the expected first-order pattern of a se- 
ries expansion in the inverse lattice volume V~ x . Reaching 
the asymptotic regime will require lattice sizes greater than 
L — 12. Our conclusions are supported by the study of many 
cumulants which all yield consistent results after proper in- 
terpretation. 



I. INTRODUCTION 

The U(l) lattice gauge model has been studied from 
the early days of numerical simulations of gauge theo- 
ries jjj . Its apparent simplicity makes it a natural choice 
for testing concepts and techniques used in simulations 
of non-abelian lattice gauge theories. However it turned 
out that the abelian character of the U(l) model brings 
specific difficulties. It can be proven rigorously under 
rather general assumptions Q that abelian lattice gauge 
theories have a phase transition at finite coupling though 
the order is not known. The determination by numerical 
simulations of the properties of the critical point in com- 
pact U(l) lattice gauge theory with Wilson action and 
periodic boundary conditions has been strongly depen- 
dent of the computing power available at the time. 

The early simulations || made runs with a few hun- 
dred iterations at each coupling constant on rather small 
(4 4 to 6 4 ) lattices. They reported a second-order phase 
transition but observed some signs of metastability with 
longer runs at isolated points. Subsequent simulations 0] 



made runs with a few thousand iterations on larger (8 4 
up to 16 4 ) lattices. It became possible to histogram the 
runs and to observe a gap in the plaquette energy. They 
concluded to the existence of a first-order critical point 
but did notice a decrease of the gap with the lattice size 
which was later confirmed Q on 9 4 and (3\/3) 4 lattices. 

But the situation became still more confusing in the 
following years. Indeed the results, always for lattices 
with periodic boundary conditions and with runs of 
O(10 4 ) length, have usually depended on which method 
was used for studying the phase transition. Most studies 
using the Monte-Carlo renormalization group transfor- 
mation |6|-|9| have claimed a second order phase tran- 
sition, but sometimes with very different critical expo- 
nents, whereas most works based on finite size scaling 
theory |]l0|dT2| have concluded to a first-order transition, 
but again with a critical exponent not consistent with 
the first-order prediction. These discrepancies led the 
subject into an expectant state. 

The interest about the nature of the phase transition 
has been revived more recently by two new kinds of re- 
sults. The role played by topological excitations, the 
monopoles, had been the object of inquiries since the 
very beginning of numerical studies of the U(l) lattice 
gauge theory |13| . The influence of topology on the crit- 
ical properties can be probed in several ways. One ap- 
proach [ l^-|l7|| is to add to the standard Wilson action a 
coupling A controlling the density of monopoles. There 
is some evidence |1| in favor of the existence of a non- 
gaussian second-order critical point in the A) plane 
when monopoles are suppressed. There is also prelimi- 
nary evidence p9| ] against universality with respect to A. 
Another approach pO| , pil| is to study the compact U(l) 
gauge theory on lattices with sphere-like topology with a 
Wilson action extended by a coupling 7 of charge 2 



S = /3^cos9 P +7^cos26 P 



(1) 



No signs of metastability are found on these lattices for 
7 < 0. A finite-size scaling analysis of the data on sphere- 
like lattices [p2| has concluded to the existence, for 7 < 0, 
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of a second-order transition with a non-gaussian contin- 
uum limit. However we must notice that the critical ex- 
ponents of these two approaches are still different. 

The latter result has spawned further investigations. 
There has been a study [E3[ of the scaling behaviour 
of gauge-ball masses and of the static potential, which 
claims to confirm the second-order nature of the transi- 
tion also on lattices with periodic boundary conditions 
at 7 = —0.2 and 7 = —0.5 with a critical exponent 
v ~ 0.36 consistent with universality and the finite-size 
scaling analysis of sphere-like lattices. The results on 
lattices with spherical topology could also lend support 
to the suspicion that the first-order signal observed on 
toroidal lattices might be a finite-size effect due to the 
topology of the lattice. This suspicion had been raised 
by the observation of the decrease of the gap mentioned 
above. However the measurements of gaps were not accu- 
rate enough and systematic enough to make a statement 
about the infinite volume limit. The particular case of 
the Wilson action was therefore re-investigated [ p4|j25| 
with increased statistics (O(10 5 ) for each coupling con- 
stant) and an emphasis on the study of gaps. There was 
definite evidence that the latent heat extrapolates to a 
non-zero value in the infinite volume limit of toroidal 
lattices. Moreover the measurement of the critical ex- 
ponent, v « 0.33, though clearly different from the first- 
order prediction, was not statistically consistent with the 
claim of universality with respect to the coupling 7. 

Therefore there is always an apparent contradiction be- 
tween the simulations on lattices with sphere-like topol- 
ogy and on lattices with periodic boundary conditions 
since one observes a gap on the latter even when 7 < 0. 
Moreover the dispersion of the predicted values for the 
critical exponent v is always as large as before, despite 
the increase in statistics since the last decade. One can 
fairly state that the confusion about the nature of the 
transition has also increased! But one can now strongly 
suspect that the origin of these paradoxical results could 
be explained by the presence of systematic corrections 
to the asymptotic finite-size scaling formulas. This will 
be our working hypothesis. There are two ways of cir- 
cumventing these corrections. The most straightforward 
approach is to simulate the largest possible lattices for a 
given computing power. However, for a transition with a 
weak first-order signal like U(l), one must thoroughly as- 
sess whether the number of tunnelling events is sufficient 
to get statistically reliable results. The other approach 
is to perform a finite-size scaling analysis of the four- 
dimensional compact U(l) gauge theory on smaller lat- 
tices but with the same quality standard as the best anal- 
yses of three-dimensional spin models. Such an analysis 
must meet two criteria. On the one hand the simulation 
must be done on many lattice sizes (and not restricted to 
3 or 4 data points as is so often the case) in order to be 
able to test the stability of the fits. On the other hand 
the statistics must be large enough to disentangle the sys- 
tematic errors from the statistical errors. Only then can 
we interpret correctly the observed deviations in a crit- 



ical exponent extracted from different observables. We 
can learn from the two and three dimensional numeri- 
cal studies of spin models how large the statistics must 
be: one needs at least O(10 6 ) configurations at each cou- 
pling constant. None of the presently published studies 
of the U(l) phase transition satisfy both constraints and 
we contend that this is the source of the contradictory 
results. 

The purpose of this work is to present a study at eight 
lattice sizes with statistics up to O(10 7 ) at each pseudo- 
critical coupling after reweighting and to clarify the na- 
ture of the U(l) phase transition. In the next section 
we shall review the part of finite-size scaling theory that 
we shall need in the sequel to interpret the data. Then 
we shall describe the details of our Monte-Carlo simu- 
lation and present the results of the measurements. In 
the following sections we shall discuss, for various cu- 
mulants, the finite-size scaling analysis of their pseudo- 
critical couplings and of their extrema. The conclusion 
will be devoted to putting all these results into a coherent 
perspective. 

II. FRAMEWORK OF THE FINITE-SIZE 
SCALING ANALYSIS 

A. Scaling ansatz 

Most of the numerical finite-size scaling analyses of 
the phase transition in the U(l) lattice gauge theory have 
been restricted to the specific-heat and the Binder cumu- 
lant. When scaling violations are important, as it turns 
out for U(l), introducing higher-order cumulants brings 
some improvement as we shall explain below. Most re- 
views on finite-size scaling theory discuss the magnetic 
cumulants only. Studying energy cumulants carries some 
specific features. Even if this is standard lore we shall 
briefly remind some useful formulas not always easy to 
find elsewhere. First, to fix our notations, we recall that 
the canonical partition function of a lattice gauge theory 
on a (i-dimensional cubic lattice of size L can be written 
as 

Z{p,L) = J n(E,L)e- /3VE dE (2) 

= e-^W) (3) 

where £l(E, L) is the microcanonical partition function, 
E is the plaquette energy and the volume, for a lattice 
gauge theory is V — \d{d — l)L d . 

Standard arguments J2(| lead for continuous phase 
transitions to the free energy density decomposition in 
a singular part and an analytic contribution: 

F(f3,L) = L- d f (x) + f ns (t,L) (4) 

where x = \t\L» is the scaling variable and t = ■£■ — 1 
is the reduced coupling. In this parameterization, the 
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hyperscaling relation a = 2 — vd is assumed with only 
one relevant scaling variable x and no dangerous irrele- 
vant variable. Moreover one assumes that there are no 
marginal variables and no logarithmic bulk singularities, 
even if we work in dimension d = 4, which could be the 
upper critical dimension of the continuum limit if there 
were a second-order phase transition. 

The function f ns represents the non-singular contri- 
bution of the background to the free energy density. 
f ns (t,L) is an analytic function in t. Its dependence 
upon L is not so clear. There is a conjecture p7| that for 
periodic systems one can take f ns (t,L) w f{t, oo). Wc 
shall assume that this conjecture holds true and that the 
L-dependence is exponentially suppressed: 



f ns (t, L) = / 00 + f 01 t + f Q2 t 2 + --- + (e- L /Z 



(5) 



However one must be aware that for systems with free 
boundary conditions one expects in general that F(f3, L) 
can have geometrical terms present, which are propor- 
tional to inverse integral powers of l/L. 

There are further analytic corrections introduced by 
non-linearities in the scaling variables away from criti- 
cality either when solving the renormalization group be- 
yond the linear approximation, t' — t + 0(t ), or for 
instance, when using t' = c //3 — 1 instead of t. Express- 
ing the non-linear scaling variable x' = t' 'L 1 ' v in terms 
of x yields corrections in powers of L~~ . We will neglect 
these corrections. 

In the renormalization group description of second- 
order phase transitions, one expects corrections to the 
asymptotic scaling function fo(x) induced by the irrel- 
evant variables. We shall keep the leading contribution 
in these corrections and parameterize it by an exponent 
lu > 0. Therefore we write 



F(J3, L) = L- d (f (x) + L-"h{x)) + f ns (t) 



(6) 



We stress that this ansatz for the free energy density is a 
natural, but completely phcnomcnological generalization 
of the asymptotic scaling formula. Moreover this ansatz 
still depends at least upon 4 parameters. It is not yet 
possible to extract that many parameters from a numer- 
ical finite-size analysis. We will have to resort to further 
approximations and limit ourselves to considering three- 
parameter fitting ansatze. 

Ansatz (||) can also describe the finite-size behavior 
at first-order transitions. Even if there is not so much 
known about this behavior from a theoretical standpoint, 



one usually expects 



that v 



d. Heuristic 



arguments based on the double-gaussian approximation 
p9|-pl| predict that the corrections should be expressible 
as a power series in V~ x . These arguments can be put 
onto a rigorous basis fi2|-|3~'|| in the special case of q-states 
Potts models for large q. 



B. Standard cumulants 

In principle we would better like to determine directly 
the finite-size scaling properties of f2(-E, L) or equiva- 
lently of the probability distribution P(E, L) of the pla- 
quette energy. In practice it is more convenient to extract 
the different moments of the plaquette distribution: 



(E n ) 

d(E n ) 
d(3 



(-iy 



1 



d n Z 



V n Z(p, L) d(3 n 
V((E) (E n )-(E n+1 )) 



(7) 
(8) 



Then, with ansatz (||) it is easy to derive the scaling 
properties of any cumulant from these moments. The 
choice of cumulants to include in the analysis is, in a 
large respect, rather arbitrary. Of course, for the sake 
of comparison with previous works, we have to study 
standard cumulants such as the specific heat: 



C v (f3,L) = -p 2 -^F(f3,L) 

= P 2 v((E 2 )-(Ef) 



and the Binder cumulant: 

We shall also introduce the second cumulant: 



U 2 (P,L) = 1 



(E) 2 



(9) 
(10) 

(11) 

(12) 



These low-order cumulants are sensitive to the analytic 
contribution to the free energy density in Eq. (^|). For 
instance ansatz (^|) implies the following for the scaling 
behavior of the specific heat: 



P 2 

C v iP,L) = -^L- 



7{'(x))-2/oo^ (13) 



But we know that \ < v < | and we can expect that 
u) i=a ^ — rf. Such an approximate equality would imply 
that the next-to-leading contribution comes both from 
the corrections to scaling and from the analytic back- 
ground. Then we can anticipate difficulties in describ- 
ing the scaling behavior of the specific heat with three- 
parameter ansatze. 

In the same way the predicted finite-size scaling be- 
havior of the second cumulant and the Binder cumulant 
is: 

U 2 (J3, L) = L- 2d+ i (u 20 {x) + u 21 (x)L- d+ i 



■OIL 



-2d+- 



,L~ 



(14) 
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U A (J3, L) = L- 2{ - d ~^ (u m (x) + U4i(x)L-( d -W 

+ 0(L- 2 ( d -h\L d -i,L-")) (15) 



The corrections to scaling are now governed by the expo- 
nent d — — . We can expect to run into problems describ- 
ing these corrections when v w 4. In the limit v = 4, 
for first-order transitions, Eqs. ( ^3| , ^4p5| ) yield the cor- 
rect leading behavior but the corrections become of order 
L- d . 

One can notice |3f| that deriving Eqs. (|l^,|l4|,|l5|) with 
respect to j3 introduces an additional factor L~ through 
the scaling variable x. Therefore studying the finite-size 
scaling behavior of the derivatives of the standard cu- 
mulants should make it easier to determine the critical 
exponent v. We have also introduced these cumulants 
into our analysis: 



dC v 
8(3 



V 



2(3 ({E 2 ) -(E) 
- )3 2 V ((£ 3 ) - 3 (E 2 ) (E) + 2 (E) 3 



(16) 



(17) 



- 2 (E A ) (E 3 ) 



(18) 



In order to make a comparison between different lat- 
tice sizes one still needs a prescription to define a common 
value of the scaling variable x. One usually uses the lo- 
cation of an extremum of a cumulant n((3, L) to define 
a universal value x* , independent of L. When there are 
several extrema it is natural to choose the extremum clos- 
est to the infinite volume critical point. Fig. [j] displays 
the plots of the cumulants C v x V^ 1 , U4 and their deriva- 



tives 



x V- 



xV 1 as a function of (3 and which 



dlh 

dp ~ ' ' 8(3 

extremum we have studied in each case. They have been 
produced by rcweighting a run with 10 6 iterations on a 
12 4 lattice at /3 = 1.01024. We recall that the critical 
coupling of U(l) lattice gauge with Wilson action and 
periodic boundary conditions is /3 c (oc) w 1.011. 



C. Derivatives of the free energy density 

We have just seen that it is not always easy to interpret 
the scaling corrections to the standard cumulants or their 
derivatives. The interpretation would be more tractable 
if we could study cumulants which are directly express- 
ible in terms of the free energy. In particular we can 



3.5C-4 
3.0e-4 
2.5o-4 
2.0e-4 
l.ao-4 
1.0e-4 



~i 1 1 1 r 

1 1 



5.06-1 
4.06-1 
3.06-1 
2.06-1 
1.0e-l 
O.OeO 
-l.Oc-1 
-2.0c-l 
-3.0c-l 
-4.06-1 

5 0c-5 ' ' ' ' ' ' ' ' -5 0c-l 

1.0094 1.0096 1.0098 1.01 1.0102 1.0104 1.0106 1.0108 1.011 



r, 



-2.0e-4 
-3.0e-4 
-4.0e-4 
-5.0e-4 
-6.0c-4 
-7.0c-4 
-8.06-4 
-9.06-4 
-1.0e-3 



~i 1 1 1 1 1 r 



-l.le-3 ' 1 1 1 ^ ' '■ — ' 1 ' -1.5C-5 

1.00941.00961.0098 1.01 1.0102 1.01041.01061.0108 1.011 

FIG. 1. The cumulants C v x V^ 1 ,Ua and their derivatives 
nk- x V ~ 1 ' ^Wjt x on a 124 lattice as a function of j3. 
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introduce the energy cumulants K n ((3,L) |3(| which are 
defined through the Taylor expansion of the free energy 
density F((3,L): 



F(p',L)=F(p,L)-J^ 



(-l) n 

L-^ K „(/?,L) (/?'-/?)" (19) 

r). 



The finite-size scaling behavior of these energy cumulants 
is much simpler since they are simply the derivatives of 
the free energy density: 



(-l) n+1 K n (p,L) 
1 



-L~ 



_ d n F((3 1 L) 

^(t\x)+L-f[ n \x))+ftH(3) (20) 



The analytic piece fi™\f3) can be neglected as soon as 
n > 3. The first three cumulants coincide with the cen- 
tral moments: 



K X = (E) 

K2 = V ((£ 2 ) - (S) 2 ) = Ma 
n z = V 2 ((E-(E)f) = Ms 



(21) 

C v /(3 2 (22) 
(23) 



The higher-order cumulants can be expressed as non- 
linear combinations of the central moments \i n — 
V n ~ 1 ((E — (E)) n ). We shall introduce in our analysis 
the cumulants K n up to sixth order: 
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FIG. 2. The cumulants « n x V~ n+1 on a 12 4 lattice as a 
function of j3. 



n± = 3V>2 (24) 

«5 =1*5- 10V/Z 2 /i3 (25) 

k 6 = - 15^/i2Ai4 - IOFM3 + 301/2 (26) 

Fig. @ displays the plots of the cumulants n n x V~ n+1 on 
a 12 lattice as a function of /3 and which extremum we 
have studied in each case. 

These energy cumulants have not been studied very of- 
ten in finite-size scaling analyses of phase transitions. In 
lower dimensions and for spin models the emphasis is of 
course on magnetic cumulants and the standard set of cu- 
mulants is natural in this context. The energy cumulants 
k„ are more suitable to the analysis of four-dimensional 
field theories. 



III. SIMULATION DETAILS 



Previous studies 1 24 2J| provide us already with a good 
knowledge of the locations of pseudo-critical couplings at 
eight lattice sizes. We have therefore decided to increase 
the statistics at several couplings inside the error bars of 
these pseudo-critical points. For each coupling we did at 
least two runs with hot and cold starts, different random 
generators with large periods, and statistics of 5 • 10 5 
configurations. For a complete overview of our statistics 
we refer to Table |. The simulated /3-values are quite 
close to each other and we can use the spectral density 
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1 x 10 6 
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1 x 10 6 
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1 x 10 6 
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1 x 10 6 


Mc 


1.01024 


4.2 x 10 6 








Mc 


1.0102 


1.4 x 10 6 



TABLE I. Number of configurations generated on different 
lattice sizes. Mc indicates multicanonical data. Simulations 
were performed in steps of 10 -4 . 



method |37| | to study the pseudo-critical points without 
noticeable extrapolation errors. 

We have of course simulated the full U(l) group in 
double precision. We used two different programs to gen- 
erate our configurations. One is an improved heatbath 
algorithm, the other one is capable of multicanonical up- 
dates (Mc) combined with over relaxation S,e3. Both 
programs were tuned to yield acceptance rates of about 
60%. The multicanonical algorithm was used in addition 
to the heatbath algorithm on 12 4 lattices. To produce 
the histogram needed for the multicanonical simulation 
we generated 1.5- 10 5 configurations using an overrelaxed 
Metropolis algorithm. For (3 — 1.01020 the multicanon- 
ically reweighted histograms of the runs at (3 = 1.01024 
were used as the input. In all Mc runs we used eight 
Metropolis updates followed by one overrelaxation step 
per link. A drastic reduction of the tunneling time is ev- 
ident in figure |[ On smaller lattice sizes the heatbath 
algorithm is sufficient to generate statistics that allow a 
precise determination of all observables. 

For safety reasons all results on our data sets were 
cross-checked with two independently developped evalu- 
ation programs. As a test of consistency we reweighted 
the multicanonical histogram to a coupling where we had 
heatbath data (see figure f|). Though we did not apply 
any smoothing algorithms to the histograms no relevant 
deviation can be observed. 

To calculate the errors of our observables we proceeded 
as follows: We first divided our runs into five bins and 
calculated the error in each run by binning. The error 
for each lattice size was then calculated by a x 2 ~ fit to a 
constant of all runs' binning results. In other words we 
did not recombine the histograms at the pseudo-critical 
points and used each run as an independent sample. Sec- 
ondly we performed a jacknife error analysis in the same 
manner but with ten bins. The larger error was taken 
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FIG. 3. The histories in multicanonical and heatbath sim- 
ulations on a 12 4 lattice and the corresponding histograms 
normalized to unit area. 
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-.565(12)E-04 


7 


1.00501(5) 


-.2404(ll)E-02 


1.00708(5) 


-.359(15)E-04 


8 


1.00728(4) 


-.1851(12)E-02 


1.00861(4) 


-.255(13)E-04 


9 


1.00851(3) 


-.1513(ll)E-02 


1.00941(3) 


-.195(12)E-04 


10 


1.00932(5) 


-.1254(14)E-02 


1.00996(4) 


-.148(9)E-04 


12 


1.010203(7) 


-.979(6)E-03 


1.010538(9) 


-.1122(12)E-04 



TABLE II. Pseudo-critical couplings and extrema of 
C v ,U2,Ua and their derivatives as a function of the lattice 
size L. 



-30 



P < £ ) 25 - 

20 - 

15 - 

10 - 

5 - 

(I - 




0.59 0.6 0.61 0.62 0.63 0.64 0.65 0.66 0.67 

E 

FIG. 4. The histograms in Mc (solid) and heatbath simu- 
lations (dashed) on a 12 4 lattice. 



into account. The jacknife error appeared to be usually 
larger, especially for the cumulants' derivatives, in heat- 
bath data at all lattice sizes except L = 12. For L = 12 
the binning error was very large in heatbath data while 
no significant difference could be seen between the two 
methods in Mc data. We interpret this phenomenon as 
a consequence of the lack of tunnelling events within a 
bin on a 12 4 lattice with a local heatbath algorithm. We 
observe just a few phase flips in 10 5 iterations. It is our 
experience that one needs at least a total of O(10 2 ) flips 
in order to control the statistical errors on the observ- 
ablcs. 



IV. MEASUREMENTS 

The measurements of the standard cumulants and their 
derivatives are displayed in Table || for all the lattice sizes 
that we have studied and the result s fo r the cumulants 
k„(3 < n < 6) are gathered in Table |lll[ 

It can be read from the tables that the statistical ac- 
curacy of the measurements of the pseudo-critical cou- 
plings and of the cumulants extrema has been improved 
by roughly one-order of magnitude with respect to pre- 
vious studies. In particular the relative accuracy of our 



G 



L 


AM 


K 3 x y- 2 




K 4 X V- 3 


4 


0.9902(2) 


-0.511(6)E-04 


0.9983(2) 


0.156(3)E-05 


5 


1.00035(6) 


-0.264(3)E-04 


1.00452(6) 


0.659(9)E-06 


6 


1.00499(4) 


-0.156(4)E-04 


1.00737(5) 


0.334(9)E-06 


7 


1.00735(5) 


-0.102(5)E-04 


1.00882(5) 


0.191(9)E-06 


8 


1.00877(4) 


-0.74(4)E-05 


1.00973(4) 


0.126(7)E-06 


9 


1.00951(3) 


-0.57(4)E-05 


1.01016(3) 


0.90(6)E-07 


10 


1.01002(6) 


-0.43(2)E-05 


1.01048(6) 


0.68(8)E-07 


12 


1.010564(9) 


-0.338(4)E-05 


1.010807(10) 


0.460(7)E-07 


L 




K S x V" 4 




k 6 x V- b 


4 


0.9851(2) 


0.385(7)E-06 


0.9907(2) 


-0.278(6)E-07 


5 


0.99827(6) 


0.131(3)E-06 


1.00110(6) 


-0.770(13)E-08 


6 


1.00381(4) 


0.55(2)E-07 


1.00543(4) 


-0.273(8)E-08 


7 


l.UUUUO^J J 




l.UU / U4i^O J 


nii 7((<w ns 

-u. ± ± I ytj J-Cj-uo 


8 


1.00830(4) 


0.162(13)E-07 


1.00895(4) 


-0.62(4)E-09 


9 


1.00919(3) 


0.11(9)E-07 


1.00964(3) 


-0.38(3)E-09 


10 


1.00979(6) 


0.71(8)E-08 


1.01011(6) 


-0.24(4)E-09 


12 


1.010446(8) 


0.452(9)E-08 


1.010611(9) 


-0.137(3)E-09 



TABLE III. Pseudo-critical couplings and extrema of k„ 
cumulants as a function of the lattice size L. 



results for the pseudo-critical couplings on the 12 4 lat- 
tice is now below 10 -5 , which is nearly comparable to 
the best results of numerical studies of lower-dimensional 
spin models. 

In the following sections we shall present a finite-size 
scaling analysis of these results. But our analysis of cor- 
rections to scaling will be conventional and far from ex- 
haustive. With our measurements given in raw form, 
without any interpretation, the reader has the possibility 
to do alternative analyses. 



V. PSEUDO-CRITICAL COUPLINGS 

The pseudo-critical couplings are expected to follow 
the asymptotic finite-size scaling formula: 



P C (L) =Pc(oo) + aL- 



(27) 



With data at eight lattice sizes it is possible to fit the 
three unknown parameters for each cumulant indepen- 
dently. With our statistics it is even possible to test the 
stability of the fits with respect to the lattice size. Ta- 
ble IV displays the results of such fits for the standard 
cumulants and their derivatives. 

The fits are clearly not stable. The statistical errors 
are small enough to show the systematic decrease of the 
critical exponent v with the lattice size L. There is also a 
slight systematic decrease of the infinite- volume limit of 
the critical coupling /J c (oo). Moreover the fits are pretty 
consistent for all cumulants. Therefore we can try a com- 
bined fit to the pseudo-critical couplings of all standard 
cumulants and their derivatives, with a common value 
of v and /3 c (oo). Such a fit is meaningful since these cu- 
mulants are defined in terms of algebraically independent 



Cumulant 


L 


x 2 


v- 1 


a 




Cv 


L>4 


2.42 


2.999(17) 


-2.091(56) 


1.01145(3) 




L > 5 


1.42 


3.037(23) 


-2.226(81) 


1.01141(3) 




L > 6 


0.731 


3.111(46) 


-2.54(21) 


1.01135(5) 




L > 7 


0.379 


3.23(11) 


-3.17(65) 


1.01128(7) 


U 2 


L > 4 


1.29 


3.090(16) 


-2.554(64) 


1.01140(2) 




L > 5 


1.04 


3.112(22) 


-2.648(92) 


1.01138(3) 




L > 6 


0.829 


3.161(44) 


-2.89(23) 


1.01134(4) 




L>7 


0.704 


3.26(11) 


-3.48(68) 


1.01128(7) 


u 4 


L > 4 


1.24 


3.116(16) 


-2.727(67) 


1.01139(2) 




L > 5 


1.08 


3.136(21) 


-2.818(97) 


1.01137(3) 




L > 6 


0.878 


3.184(43) 


-3.07(24) 


1.01133(4) 




L > 7 


0.816 


3.28(11) 


-3.66(71) 


1.01127(7) 


96V 

9/3 


L > 4 


1.43 


3.012(26) 


-1.398(58) 


1.01136(3) 


L > 5 


1.01 


3.057(37) 


-1.509(91) 


1.01133(3) 




L > 6 


0.715 


3.140(71) 


-1.75(22) 


1.01128(5) 




L > 7 


0.730 


3.26(17) 


-2.22(72) 


1.01123(7) 


t)U 2 

9/3 


L > 4 


1.12 


3.096(24) 


-1.705(66) 


1.01133(2) 


L > 5 


0.836 


3.133(35) 


-1.81(11) 


1.01131(3) 




L > 6 


0.704 


3.196(67) 


-2.03(24) 


1.01128(4) 




T > 7 


u.oyo 


q OQ/1 R\ 

O.ZO ^ 




i ni i r >A(7\ 


as 


L>4 


0.724 


3.142(24) 


-1.904(71) 


1.01131(3) 


L > 5 


0.687 


3.165(34) 


-1.98(11) 


1.01130(3) 




L > 6 


0.613 


3.218(66) 


-2.18(26) 


1.01127(4) 




L > 7 


0.732 


3.30(16) 


-2.57(77) 


1.01123(7) 



TABLE IV. Independent second-order finite-size scaling 
fits to the pseudo-critical couplings of standard cumulants 
and their derivatives. 



combinations of the moments of the plaquette energy dis- 
tribution. We can also test the stability of this combined 
fit with respect to the lattice size. The parameters of 
these fits are given in Table 0. Fig. || displays the result 
of the combined fit to the L > 6 data points. Errors on 
the data points are much smaller than the marker sizes. 

We can repeat the same kind of analysis for the en- 
ergy cumulants K n . The independent fits are shown in 
Table VI and the combined fits in Table VII. The re- 



sult of the combined fit to the L > 6 data points is also 
displayed in Fig. [| 

The critical values of v and /3 c (oo) extracted from the 
cumulants n n are very similar to those extracted from 
the standard cumulants or their derivatives. Studying 



L 


x 2 


/3 c (oo) 


v- 1 


L>A 


5.86 


1.011348(8) 


3.1059(72) 


L > 5 


4.09 


1.011317(10) 


3.144(10) 


L > 6 


1.81 


1.011269(14) 


3.222(19) 


L > 7 


0.692 


1.011224(19) 


3.329(35) 


L > 8 


0.799 


1.011235(29) 


3.304(62) 


L > 9 


0.872E-01 


1.011206(43) 


3.4249(69) 



TABLE V. Combined asymptotic finite-size scaling fits to 
the pseudo-critical couplings of the cumulants C v , U2, U4, 
and 
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FIG. 5. Combined fit to the L > 6 data points of the 
pseudo-critical couplings of the standard cumulants and their 
derivatives. 



Cumulant 


L 


x 2 


t/" 1 


a 


/3 c (oo) 


^3 


L > 4 
L > 5 
L > 6 
L> 7 


1.72 
1.05 
0.720 
0.711 


3.005(27) 
3.060(38) 
3.146(72) 
3.28(17) 


-1.381(58) 
-1.515(92) 
-1.77(23) 
-2.26(74) 


1.01136(3) 
1.01132(3) 
1.01128(5) 
1.01122(7) 




L > 4 
L > 5 
L > 6 
L > 7 


1.18 
0.872 
0.724 
0.812 


3.011(44) 
3.078(62) 
3.20(13) 
3.38(28) 


-0.858(60) 
-0.959(97) 
-1.19(27) 
-1.69(93) 


1.01129(3) 
1.01127(3) 
1.01123(5) 
1.01119(7) 


K5 


L > 4 
L > 5 
L > 6 
L > 7 


1.02 
1.08 
0.684 
0.721 


3.083(22) 
3.063(31) 
3.139(60) 
3.24(15) 


-1.876(65) 
-1.815(92) 
-2.08(22) 
-2.51(68) 


1.01133(3) 
1.01134(3) 
1.01130(5) 
1.01125(8) 




L > 4 
L > 5 
L > 6 
L > 7 


0.843 
0.887 
0.485 
0.489 


3.103(28) 
3.079(40) 
3.174(77) 
3.29(19) 


-1.512(67) 
-1.454(95) 
-1.72(24) 
-2.137(64) 


1.01129(3) 
1.01130(3) 
1.01126(5) 
1.01122(7) 



TABLE VI. Independent second-order finite-size scaling 
fits to the pseudo-critical couplings of the k„ cumulants. 



L 


x 2 




v- 1 


L > 4 


2.41 


1.011299(12) 


3.084(14) 


L > 5 


1.59 


1.011288(13) 


3.105(18) 


L > 6 


0.718 


1.011244(17) 


3.202(31) 


L > 7 


0.442 


1.011206(25) 


3.319(63) 


L > 8 


0.564 


1.011220(37) 


3.27(11) 


L > 9 


0.856E-02 


1.011191(53) 


3.429(14) 



1.015 



1.01 



1.005 



0.995 



0.99 



0.985 



0.98 



0.975 



0.97 
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TABLE VII. Combined asymptotic finite-size scaling fits 
to the pseudo-critical couplings of the cumulants K3, K4, K5 
and kg. 



FIG. 6. Combined fit to the L > 6 data points of the 
pseudo-critical couplings of the cumulants k„. 



the pseudo-critical couplings of the cumulants K n does 
not bring any significant improvement over the standard 
cumulants. The pseudo-critical couplings are not very 
sensitive to the corrections to asymptotic scaling. Even 
if the finite-size scaling violations are clearly visible on 
Figs. [5] and |] at the small lattice sizes L = 4 or 5, it 
would not be possible to introduce additional parameters 
in the fits. 



VI. CUMULANTS' EXTREMA 

In this section we will present some selected analyses 
of the finite size scaling behaviour of the cumulants we 
calculated. As the main result we can state that the 
finite-size scaling behavior of all cumulants is consistent 
with a first-order transition. However we observe a differ- 
ence in the corrections to scaling of the specific heat and 
the Binder cumulant on one hand and their derivatives 
and the energy cumulants on the other hand. Up to the 
lattice sizes we calculated the first do not follow yet the 
expected pattern in a first order phase transition while 
the other show a clear first order behaviour. To demon- 
strate this we will present selected fits to our data. In 
all pictures in the following subsections the solid curves 
denote the fit to the ansatz which yields the best results. 



A. Asymptotic scaling 

It is certainly not possible to describe the data without 
taking into account the corrections to scaling. However 



8 



L 




Va 


90 


au A 

8/3 


9 
8 
7 
6 
5 
4 


3.299(17) 
3.177(28) 
3.121(16) 
3.030(12) 
2.983(7) 
2.968(5) 


3.250(17) 
3.132(27) 
3.070(17) 
2.942(12) 
2.856(7) 
2.788(5) 


3.467(56) 
3.19(11) 
3.23(11) 
3.123(64) 
3.050(36) 
3.016(18) 


3.402(63) 
3.19(12) 
3.183(97) 
3.064(59) 
2.962(32) 
2.906(16) 


L 






Kb 


k 6 


g 
8 
7 
6 
5 
4 


3.479(63) 
3.19(11) 
3.22(12) 
3.125(66) 
3.048(38) 
3.019(20) 


3.422(57) 
3.30(13) 
3.247(80) 
3.139(51) 
3.075(30) 
3.046(20) 


3.405(55) 
3.27(12) 
3.282(90) 
3.135(59) 
3.055(36) 
3.038(18) 


3.416(47) 
3.30(11) 
3.246(62) 
3.126(39) 
3.061(22) 
3.045(15) 



TABLE VIII. Critical exponent v~ extracted from lattice 
sizes L, L + 1 and L + 2 for each cumulant extremum with 
the ansatz FSS2a. 



it is instructive to define, for any cumulant k, an "effec- 
tive critical exponent" v K {L) by fitting the data at scales 
L, L + 1 and L + 2 to the asymptotic form of the scaling 
ansatze of sect. Q: 



FSS2a 



k{L) = bh T 



(28) 



We shall refer to this ansatz as FSS2a. The relation be- 
tween the critical exponents r and v depends of course 
up on th e cumulant. The results are displayed in Ta- 
ble Ivm. 



We observe again in Table VIII a systematic decrease 
of the critical exponent v with the lattice size L. The 
effective critical exponent v{L) is pretty consistent at 
each scale L across the cumulants n n and with the val- 
ues found in the combined fits of the pseudo-critical cou- 
plings. There is however some dispersion among the stan- 
dard cumulants and their derivatives not observed in the 
analysis of pseudo-critical couplings. In fact the x 2 01 
these asymptotic fits is quite high even with 2 parameters 
and only 3 data points. We see here a first manifestation 
of the expected improved scaling behavior of the energy 
cumulants k„. 



B. Corrections to scaling 

The asymptotic scaling fits for the cumulants' extrema 
as well as for the pseudo-critical couplings hint at a first- 
order transition. Our working hypothesis will be to de- 
scribe the data with scaling corrections to a first-order 
transition. These corrections are expected to be express- 
ible as a series expansion in the inverse volume V . 
Since we limit ourselves to three-parameter fits, we shall 
introduce the two fitting ansatze: 



FSSla 
FSSlb 



k(L) x V~ k = 

k(l) x y- k = 



-bL- d 
-bL- d 



(29) 

-cL- 2d (30) 




FIG. 7. Fits to the specific heat and its derivative. 



The volume normalization factor is chosen such that the 
fitting ansatz makes sense. With the definitions of sect. 
we have k = for U 2 and f/ 4 , k = 1 for C v , , 

k = 2 for and k = n — 1 for K n . 

We shall need another, more phenomenological fitting 
ansatz: 



FSS2b : 



k(L) x V~ k =a + bL- 



(31) 



The volume normalization factor V~ k is by definition the 
same as for the first-order fitting ansatze. Therefore the 
exponent tu parameterizes the corrections to an asymp- 
totic first-order behaviour. Introducing such an effective 
exponent allows a more flexible description of the correc- 
tions. One should find that to — ► d when we reach the 
asymptotic regime. 

The best FSSlb fit to the specific heat is represented 
by the dashed line in Fig. ^. Obviously it cannot describe 
the smallest lattice sizes. On the other hand the fitting 
ansatz FSS2b describes our data on the specific heat 
perfectly. Our data on dC v /d(3 have the largest errors 
of all the results we produced. Thus the quality of the 
best fits with FSSlb (solid line) or FSSla (dashed line) 
is not good. 

The ansatz FSS2b is also able to reproduce the data 
on the Binder cumulant quite accurately (solid line in 
Fig. ||). We observe again that ansatz FSSlb (dashed 
line) fails to reproduce the smallest lattice sizes. The 
description of the derivative dXJ^jdfi by FSSlb (solid 
line) is much better and FSSla gives already a good fit 
for all but the smallest lattice sizes (dashed line). The 
results for U2 and its derivative are quite similar to U4 
and its derivative and will not be reproduced here. 

In the same way we find that the first-order ansatze 
can fit the data on the energy cumulants K n . In Fig. ^ 
all solid lines represent the best fits with ansatz FSSlb 
whereas all dashed lines represent the best fits with 
ansatz FSSla. The cumulant K4 is a noticeable ex- 
ception. The asymptotic first-order ansatz FSSla is al- 
ready able to describe all data points and is nearly in- 
distinguishable from ansatz FSSlb. We have another, 
and more vivid manifestation of the improved scaling be- 
haviour of the energy cumulants K n . 
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FIG. 8. Fits to the Binder Cumulant and its derivative. 
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FIG. 9. Fits to the energy cumulants k„. 



Cum. 


Fit 


L 


x 2 


a 


b 


C v 


FSS2b 


L > 5 


0.52 


.147(6)E-03 


.580(23)E-01 


a/3 


FSSlb 


L > 5 


0.58 


.149(4)E-04 


.121(5) 


Ui 


FSS2b 


L > 4 


0.85 


-.545(ll)F-03 


-.352(7) 


au 4 

8/3 


FSSlb 


L > 4 


1.79 


-.823(14)E-05 


-.632(11)E-01 




FSSlb 


L > 5 


0.50 


-.244(6)E-05 


-.196(9)E-01 


K4 


FSSla 


L > 4 


0.94E-01 


-.270(8)E-07 


-.395(5)E-03 


«5 


FSSlb 


L > 4 


1.26 


.143(15)E-08 


.631(22)E-04 


K6 


FSSlb 


L > 5 


0.43 


-.334(71)E-10 


-.206(13)E-05 



TABLE IX. Best fits to the cumulants' extrema. 



Table IX gives a quantitative content to the previous 
qualitative observations. This table contains the param- 
eters a and b of the best fits displayed in the figures. 
We always choose the smallest lattice size with a x 2 P er 
degree of freedom less than 2. 

The main feature is the smallness of the ratios f . In 

b 

the first-order fits, the leading contribution a is generally 
smaller than the correction term bL~ d up to the lattice 
size L = 12. It is no surprise that the first-order nature 
of the transition is so difficult to observe. 

We can get order-of-magnitude estimates of the param- 
eters of a first-order transition by comparing the values of 
a with the predictions of the (very crude) double-gaussian 
approximation p9f . The height of the maximum of the 
specific heat should increase linearly with L d as: 



O(l) 



(32) 



where /3 C is the infinite-volume critical coupling and e 
e<z is the latent heat. We find 



e Q - e d « 0.029 



(33) 



The minimum of the Binder cumulant is predicted to be 
pOpl: 



{el -elf 
12(e e rf ) 2 



+ o(v- 1 ) 



(34) 



We get another estimate of the latent heat which agrees 
pretty well with (fj3|): 

e - e d « 0.026 (35) 



VII. CONCLUSION 

Our study shows that finite-size scaling violations are 
indeed present in each observable for lattice sizes 4 < 
L < 12. The observed scaling violations are consistent 
for all observables: the critical exponent v systematically 
decreases with the lattice size L. Asymptotic finite- 
size scaling of all cumulants and of their pseudo-critical 
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couplings yields consistently v w 0.29 for L > 9 which 
points towards a first-order transition. 

The scaling violations in the pseudo-critical couplings 
and the cumulant values decrease slowly with the lat- 
tice size. This slow variation, which is hard to unravel, 
can explain the claims for a second-order transition with 
v w 0.33. However the scaling violations seem to de- 
crease more rapidly for the derivatives of the standard 
cumulants and for the cumulants k„. The finite-size size 
behavior of all these cumulants can be completely de- 
scribed in the range 4 < L < 12 by volume correction 
terms of order V~ l and V~ 2 . The cumulant K4 is even 
completely described by the asymptotic first-order for- 
mula. 

There seems to be a correlation between the amount 
of scaling violations in each cumulant extremum and the 
location of its pseudo-critical-coupling. The closer the 
pseudo-critical coupling is, at fixed lattice size, from the 
infinite- volume critical point, the better the description 
of the extremum by a first-order transition. All the cu- 
mulants we have studied follow this rule. Then we can 
make two observations. 

On the one hand, if we had chosen to study the ex- 
tremum in each cumulant which is farthest from /3 c (oo) 
we would certainly have found a behavior inconsistent 
with a first-order transition and concluded to a second- 
order transition with a critical exponent v greater than 
0.33. We suggest that this is the origin of many claims for 
a large critical exponent v. Most of them are not based 
on finite-size scaling analyses of the standard cumulants, 
which all should give v < 0.33, but come from various 
analyses done at couplings f3 smaller than f3c v {L), the 
pseudo-critical coupling of the specific heat in the corre- 
sponding lattices. 

On the other hand we could think of studying still 
higher-order k„ cumulants to get closer to /3 c (oo). How- 
ever we would face two difficulties. First the statistical 
noise increases with n. Secondly the extrapolation from 
the couplings where we generated our configurations gets 
large. The solution might be to generate the configu- 
rations for all lattice sizes at the infinite-volume critical 
point and use the method developed [|2|3^] for the study 
of first-order transitions in Potts models. 

However it is not clear how the proofs can be extended 
to the case of the U(l) phase transition. Potts models 
have a discrete symmetry and a local order parameter 
whereas the U(l) lattice gauge theory has a continuous 
local symmetry and a non-local order parameter since the 
proof of existence of the phase transition uses the Wilson 
criterion |^7| . In Potts models the physical correlation 
lengths in both pure phases stay finite in the infinite vol- 
ume limit whereas there exist massless photons in the 
ordered phase of U(l). Clearly the U(l) lattice gauge 
theory deserves as much numerical study as the lower- 
dimensional spin models and theoretical understanding 
of the scaling violations will be required before a defini- 
tive conclusion on the nature of its phase transition. 



ACKNOWLEDGMENTS 

One of us (OR.) is grateful to the Freie Univer- 
sitat Berlin for its hospitality during a visit when this 
manuscript was completed and acknowledges a financial 
support by the Graduiertenkolleg "Strukturuntersuchun- 
gen, Prazisionstests und Erweiterungen des Standard- 
modells der Elementarteilchenphysik" . 

B. K. wants to thank P. E. L. Rakow from DESY-IfH 
Zcuthcn for his support during the preparation of the 
diploma thesis. He initialized the investigation of the 
latent heat. 



[1] M. Creutz, L. Jacobs and C. Rebbi, Phys. Rev. D20, 
(1979) 1915. 

[2] J. Frohlich and T. Spencer, Commun. Math. Phys. 83 
(1982) 411. 

[3] B. Lautmp and M. Nauenberg, Phys. Lett. B95 (1980) 
63. 

[4] J. Jersak, T. Neuhaus and P.M. Zerwas, Phys. Lett. 

133B (1983) 103. 
[5] R. Gupta, M.A. Novotny and R. Cordery, Phys. Lett. 

B172 (1986) 86. 
[6] A. N. Burkitt, Nucl. Phys. B270 [FS16] (1986) 575. 
[7] C. B Lang, Nucl. Phys. B280 [FS18] (1987) 255. 
[8] C.B. Lang and C. Rebbi, Phys. Rev. D35 (1987) 2510. 
[9] K. Decker, A. Hasenfratz and P. Hasenfratz, Nucl. Phys. 
B295 [FS21] (1988) 21. 
[10] H.G. Evertz, J. Jersak, T. Neuhaus and P.M. Zerwas, 

Nucl. Phys. B251 [FS13] (1985) 279. 
[11] J. Jersak, T. Neuhaus and P.M. Zerwas, Nucl. Phys. 

B251 [FS13] (1985) 299. 
[12] G. Bhanot, T. Lippert, K. Schilling and P. Ueberholz, 

Nucl. Phys. B378, (1992) 633. 
[13] T. DeGrand and D. Toussaint, Phys. Rev. D22 (1980) 
2478. 

[14] J.S. Barber, Phys. Lett. B147 (1984) 330. 

[15] J.S. Barber, R.E. Shrock and R. Schrader, Phys. Lett. 

B152 (1985) 221. 
[16] J.S. Barber and R.E. Shrock, Nucl. Phys. B257 (1985) 

515. 

[17] W. Kerler, C. Rebbi and A. Weber, Phys. Rev. D50 
(1994) 6984; Nucl. Phys. B450 (1995) 452; Phys. Lett. 
B348 (1995) 565; Phys. Lett. B380 (1996) 346. 

[18] W. Kerler, C. Rebbi and A. Weber, Phys. Lett. B392 
(1997) 438. 

[19] G. Damm a nd W. Kerler, HUB-EP-97/63 (1997), \mp- 
lat/970906l| . 

[20] C.B Lang and T. Neuhaus, Nucl. Phys. B431, (1994) 
119. 

[21] J. Jersak, C.B Lang and T. Neuhaus, Phys. Rev. Lett. 

77 (1996) 1933. 
[22] J. Jersak, C.B Lang and T. Neuhaus, Phys. Rev. D54 

(1996) 6909. 



11 



[23] J. Cox, W. Franzki, J. Jersak, C.B. Lang, T. Neuhaus, 

A. Seyfried, P. W. Stephenson, Nucl. Phys. B499 (1997) 
371. 

[24] B. Klaus, Diploma thesis (FU Berlin, 1995). 

[25] C. Roiesnel, Phys. Lett. B405 (1997) 126. 

[26] M.E. Fisher, in Critical phenomena, Proc. 51st En- 
rico Fermi Summer School, Varena, ed. M.S. Green 
(Acadamic Press, NY, 1972). 

[27] V. Privman, P.C. Hohenberg and A. Aharony, Phase 
Transitions and Critical Phenomena (ed. C. Domb and 
J.L. Lebowitz), Vol. 14, p. 1, Academic Press, (1991) 
London. 

[28] M.E.Fischer and N. Berker, Phys. Rev. B26 (1982) 2507. 
[29] M.S. Challa, D.P. Landau, K. Binder, Phys. Rev. B34 
(1986) 1841; 

[30] A. Billoire, R. Lacaze, A. Morel, S. Gupta, A. Irback and 

B. Petersson, Phys. Rev. 42 (1990) 6743. 

[31] J. Lee and J.M. Kosterlitz, Phys. Rev. B43 (1991) 3625. 
[32] C. Borgs and R. Kotecky, J. Stat. Phys. 61 (1990) 79; 
[33] C. Borgs, R. Kotecky, and S. Miracle-Sole, J. Stat. Phys. 
62 (1991) 529. 

[34] A. Billoire, T. Neuhaus and B. Berg, Nucl. Phys. B413 

(1994) 795. 
[35] K. Binder, Z. Phys. B 43 (1981) 119. 
[36] W. Janke and S. Kappler, J. Phys. I (France) 7 (1997) 

663. 

[37] A.M. Ferrenberg and R.H. Swendsen, Phys. Rev. Lett. 

61, (1988) 2635; 63 (1989) 1658; 63 (1989) 1195. 
[38] B. Berg and T.Neuhaus, Phys. Lett. B267 (1991) 249; 

Phys. Rev. Lett. 68 (1992) 9. 



12 



